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ABSTRACT 

We analytically investigate the condition for a particle accelerator to be active in the outer mag- 
netosphere of a rotation-powered pulsar. Within the accelerator (or the gap), magnetic-field-aligned 
electric field accelerates electrons and positrons, which emit copious gamma-rays via curvature pro- 
cess. If one of the gamma-rays emitted by a single pair materializes as a new pair on average, the 
gap is self-sustained. However, if the neutron-star spin-down rate decreases below a certain limit, 
the gap becomes no longer self-sustained and the gamma-ray emission ceases. We explicitly compute 
the multiplicity of cascading pairs and find that the obtained limit corresponds to a modification of 
previously derived outer-gap death line. In addition to this traditional death line, we find another 
death line, which becomes important for millisecond pulsars, by separately considering the threshold 
of photon-photon pair production. Combining these traditional and new death lines, we give predic- 
tions on the detectability of gamma-ray pulsars with Fermi and AGILE. An implication on the X-ray 
observations of heated polar-cap emission is also discussed. 

Subject headings: gamma rays: stars — magnetic fields — methods: analytical — stars: neutron 



1. INTRODUCTION 

The Fermi Large Area Telescope (LAT) provides a 
wealth of new data on isolated, rotation-powered pulsars, 
increasing the number of detected 7-ray pulsars from 
seven to more than sixty (e.g., Abdo et al. 2010a). The 
A GILE has also reported the detection of about twenty 
7-ray pulsars (e.g., Pellizzoni et al. 2009). Since inter- 
preting 7-rays should be less ambiguous compared with 
reprocessed, non-thermal X-rays, The 7-ray pulsations 
observed from these objects are particularly important 
as a direct signature of basic non-thermal processes in 
pulsar magnetospheres, and potentially should help to 
discriminate among different emission models. 

In a pulsar magnetosphere (fig. [1]), there is a sur- 
face called the 'null surface', on which the Goldreich- 
Julian charge density (Goldreich & Julian 1969) pgj = 
B ■ Sl/(2nc) changes sign, where B denotes the local 
magnetic field, f2 the NS angular-velocity vector, and c 
the speed of light. There is another characteristic surface 
called the 'light cylinder' beyond which the co-rotational 
velocity exceeds c. The distance of the light cylinder 
from the rotation axis is called the 'light-cylinder radius', 
?37lc = c/f2, where £1 = \Sl\ denotes the NS rotational 
angular frequency. On each magnetic azimuthal angle 
(measured around the magnetic axis), there is a magnetic 
field line that crosses the light cylinder tangentially; they 
are called the 'the last-open field lines'. 

In the lower colatitudes than the last-open field lines, 
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field lines closes inside the light cylinder; thus, a high 
plasma density is expected by the trapping of relativis- 
tic charges due to magnetic mirrors, preventing the oc- 
currence of a strong magnetic-field-aligned electric field, 
E», in this closed zone. In the higher colatitudes, on the 
other hand, field lines open to a distant region beyond the 
light cylinder, allowing a pulsar wind to flow along them, 
thereby resulting in a high vacuum with unscreened E» 
in some regions in this open-field-line zone. 

On these grounds, in all the pulsar emission models, 
particle acceleration is assumed to realize in the open 
zone. In polar-cap (PC) models, emission takes place 
within a few neutron star (NS) radii above a PC surface 
(Arons & Scharlemann 1979; Daugherty & Harding 1982, 
1996). However, the observed cutoff energies, typically 
between 1 GeV and 4 GeV (Abdo et al. 2010a), along 
with the 25 GeV detection of pulsed signals from the 
Crab pulsar (Aliu at al. 2008), suggests no attenuation 
from one-photon absorption by a strong magnetic field 
(Baring 2004) near the NS, indicating the minimum alti- 
tudes of a few stellar radii above the PC surface. Thus, 
the possibility of a high-altitude emission gathered at- 
tention. 

Recent pulsar high-energy emission models adopt, 
therefore, higher-altitude geometries. There are two 
main scenarios in this approach: the outer-gap (OG) 
model (Cheng et al. 1986a, b; Romani 1996; Cheng et 
al. 2000; Hirotani 2008; Takata et al. 2008; Romani 
& Watters 2010), and the pair-starved polar-cap model 
(Frackowiak & Rudak 2005; Harding et al. 2005; Venter 
et al. 2009). 

In the present paper, we focus on the OG model, 
in which pairs are created in the outer magnetosphere 
mainly by photon-photon (7-7) pair production. The 
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Fig. 1. — Side view of magnetic field lines on the meridional 
plane in which both the magnetic and rotational axes reside. 

produced pairs polarize owing to the magnetic-field- 
aligned electric field, E\\ . If the rotation and mag- 
netic axes reside in the same hemisphere, a positive E\\ 
is exerted to accelerate e + 's (or e~'s) outward (or in- 
wards), increasing the real charge density outwards. The 
inward-migrating, relativistic e~'s radiate curvature 7- 
rays, some of which (nearly head-on) collide with the X- 
rays emitted from the NS surface to materialize as pairs. 

The purpose here is to examine the detectability of the 
7-ray photons from pulsars with OG accelerators. On 
the neutron-star period (P) versus period-derivative (P) 
plane, we find two kinds of death lines below which the 
pulsar OGs become inactive. In § [2] we derive the first 
kind death line, which turns out to be a modification of 
the previous death lines (Zhang et al. 2004). We then 
derive the second kind, new death line in § [3j by sepa- 
rately considering the threshold of 7-7 pair production. 
We finally combine these two death lines and discuss the 
detectability of pulsars with recent 7-ray telescopes. 

2. DEATH LINE DERIVED BY MULTIPLICITY 
CONSTRAINT 

According to the analytical consideration in classic OG 
model (Cheng et al. 1986a, b; Romani 1996; Cheng et 
al. 2000), and to the numerical solution of the OG (Hi- 
rotani 2006a, b), an OG is found to exist essentially be- 
tween the null surface and the light cylinder (fig. [1]). 
Since most of the pairs are produced in the inner-most 
region of the OG, the (outward- moving) e + 's propagate 
longer distance than the (inward- moving) e~'s. As a re- 
sult, the outward 7-ray flux becomes greater than the 
inward one. This result is particularly important when 
we predict the high-energy emission properties from an 
OG. Nevertheless, to investigate the condition for an OG 
to be self-sustained, it is crucial to consider the collisions 
between the inward 7-rays, which are emitted by e _ 's, 
and the outward-propagating X-rays, which are emitted 
from the heated PC surface. Thus, in this paper, we 
concentrate on the condition for the inward-propagating 
7-rays to materialize as pairs in the inner magnetosphere. 

For young pulsars like Crab or Vela, 7-7 pair produc- 
tion is sustained by the strong thermal emission from the 
cooling NS surface, and the gap trans-field thickness is 



kept thin. However, as the NS ages, the heated PC emis- 
sion due to the bombardment of in- falling e~'s dominates 
the cooling NS emission and controls the gap trans-field 
thickness. Therefore, when we investigate the death lines 
below which an OG accelerator ceases to emit 7-rays, it 
is appropriate to consider the heated PC as the source 
of X-rays for 7-7 pair production. For simplicity, we as- 
sume that the heated PC has a uniform temperature kT 
with area 7rPp C . When the OG activity is marginally 
sustained, pair production takes place mainly near the 
null surface. 

2.1. Optical depth of photon-photon pair production 

At an altitude r — ^> P pc , the photon energy flux is 
given by 

F = fJ-^-) cos 8, (1) 

\r — r*J 

where the flux measured at the surface becomes 
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P* = 1.0 x 10 



20 



/ kT 



V 100 eV 



ergs 'cm 2 ; (2) 



r* denotes the NS radius, O the colatitude of the point 
from the magnetic axis. In this paper, we consider that 
the magnetic axis is moderately or highly inclined with 
respect to the rotation axis. In this case, the PC surface 
is almost face-on from the gap inner boundary, which is 
located well inside the light cylinder. Thus, we adopt 
cos O pa 1 in what follows. 

Evaluating the photon number flux at typical black- 
body photon energy, hv x rj 2.8fcT, we obtain the mean- 
free path for an inward-propagating 7-ray to materialize 
in a collision with the surface X-ray as 



1 C ftf» 



9.9 x 10 5 



kT 
100 eV 



Po 



(3) 



where the 7-7 pair-production total cross section is eval- 
uated by (Tp ~ 0.2(Tt = 1.3 x 10~ 25 cm 2 , because the 
collisions take place mostly head-on. 

Inward-propagating 7-rays are typically emitted from 
the central colatitude in the gap. Therefore, before cross- 
ing the critical field lines on which the null surface inter- 
sects the light cylinder (i.e., before escaping from the 
gap), they propagate the distance 
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[null 

V2 



(sin 6 nu \\ — tan a cos nu ii) 



(4) 



where r nu n denotes the distance of the null surface on the 
last-open field line from the NS center, a the angle be- 
tween rotation and magnetic axes. If a = 75° , 9 nu \\ « 85° 
gives I w 0.67r nu n, whereas if a = 45°, 8 nu n « 74° gives 
I « 0.47r nu n, for example. Note that the gap is sup- 
posed to exist only below the critical field lines, because 
such solutions are being obtained by three-dimensional 
OG simulations. Because of this geometry, the obtained 
death lines give tighter constraints than previous ones, 
which are obtained by assuming that the OG occupies 
the entire open-field-line regions between the null sur- 
face and the light cylinder (i.e., by imposing / < 1 as in 
Zhang et al. 2004). 
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Combining equations ([3]) and Q , we can compute the 
optical depth r = l/X p for 7-7 pair production, 



1.0 



/ kT 



VlOOeV 



(n-iy 



?"null 



(5) 



where rj = r/r*, r» g = ?~/(10 6 cm), P = 2ir/fl is the NS 
spin period; r rj r nu n is used. 

2.2. Acceleration electric field 

For the pulsars near the death line (to be derived be- 
low), the potential drop in the OG becomes a good frac- 
tion of the electromotive force, K ~ -B*r 3 /(2n7£ c ), ex- 
erted on the PC surface. Thus, we can reasonably eval- 
uate the maximally possible Eu by 



E 



0.5K 



(G) 



where fi denotes the NS dipole moment, and _B* the mag- 
netic field strength at the NS surface. Since the poten- 
tial drop occurs only in a portion of the open zone, we 
multiply the factor 0.5 in the numerator. By this accel- 
eration electric field, the inward motion of e^'s attain 
a force balance between the electrostatic and curvature- 
radiation forces. Thus, the typical Lorentz factors of the 
e _ 's become 

■ml ^ /4 
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2e 



-E 



(7) 



where e refers to the magnitude of the charge on the 
electron, R c the curvature radius of the magnetic field 
line. 

2.3. Evaluation of curvature radius 

To evaluate R c , we we consider the following two rep- 
resentative cases. 

Case 1: The null surface intersects the last-open field 
line in the inner-most magnetosphere; that is, r nu n < 
0.05tulc, In this case, the curvature radius can be ap- 
proximated as 



Rc 



(8) 



on the last-open field line. If a > 70°, this case becomes 
applicable along the magnetic field lines that crosses the 
NS surface with magnetic azimuthal angle in the range 
—60° < <p* < 60°, where increases counter-clockwise 
and ip* = 7r points the rotation axis from the magnetic 
axis. To derive the OG death line, it is sufficient to 
consider only this range, |y*| < 60°, because the OG is 
active only in \ip*\ < 60°. 

Case 2: The null surface intersects the last-open field 
line in the outer magnetosphere; that is, r nu n > 0.07tulo 
In this case, we have 

R c w 0.5tu LC . (9) 
If a < 62°, this case becomes applicable in \(p*\ < 60°. 

2.4. Typical gamma-ray energies 
Substituting equation (|6|) into 0, we obtain 



7 «3.3xl0V/4p-i/2 M3o i/4 r ^i/4 j 



(10) 



for case 1, and 



7 ; 



1.7 x 10 7 P" 1/4 



fJ-30 



1/4 



(11) 



for case 2, where = n/ (10 30 G cm 3 ). 

The flux density of the curvature spectrum peaks at the 
frequency 0.29^ c , where the critical frequency is defined 
by 

-ft <-> 

However, for the pulsars near the death line, 7-7 pair 
production is maintained by the collisions between the 
7-ray photons with energy hv 1 3> 0.29hv c and the ther- 
mal photons near the blackbody peak hv s ~ 2.82fcT, or 
between the 7-rays with hv 7 ~ 0.29/i^ c and the thermal 
photons in the Wien regime hv s ^> 2.82/cT, or both (i.e., 
hv^ > 0.29hi/ c and hv s > 2.82&T). To take account of 
such effects, we can practically assume that pair produc- 
tion takes place between the 7-rays with hv~ ~ hv c and 
the surface X-rays with hv s ~ 3fcT; this treatment in- 
curs relatively small errors from the exact, complicated 
numerical computations. 

Substituting equations (fT0|) or (fTTj) into (fl2")) , we obtain 
the typical 7-rays energy, 



hv c = 11.8 r) 1/4 P~ 



W /4 r*, 6 1/4 MeV, 



hv c = 60P" 7/ V30 3/4 MeV 
for cases 1 and 2, respectively. 



(13) 
(14) 



2.5. Number of gamma-rays emitted by each electron 

Let us now evaluate number of photons emitted by a 
single inward-accelerated e~ near the null surface, the 
inner-most region of the gap. If R c were constant in 
the magnetosphere, the photons emitted between the al- 
titudes r nu n and (r nu n + I) illuminate the null surface, 
where I is given by equation Therefore, we could 
evaluate the number of 7-rays emitted by a single e~ by 

N 1 » — J-. (15) 



It follows that a single e~ cascade into N 7 t secondary 
pairs within the gap. 
Substituting equations fp~3|) or (fl~4|) into (fl~5|) . we obtain 



7V 7 «78 r, 3 / 4 p-V 3 o 1/4 n, 6 3 / 4 



7V 7 « 10.8 77^~ 5/4 M30 1/4 r*, 6 



for cases 1 and 2, respectively. 



0.67r nu n ' 
I 



0.47r n 



(16) 
(17) 



2.6. Number of secondary pairs cascaded from a 
primary electron 

Combining equations ([5]) and (|16p , we obtain the num- 
ber of pairs that a single e~ cascades within the gap, 



iV 7 r«52 



r) 



7/4 



(V I) 2 



100 cV 



xrV 3 » 1/4 r,; /4 



0.67r nuU 



(18) 
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for case 1. In the same manner, combining equations ([3]) 
and (HTJ), we obtain 



iV T T~ 5.0 



T 



(r, l)^ 

x P- 5 ' A 



R 



pc 



kT 
100 eV 
I 

0.47r null 



(19) 



for case 2. We have to express (R pc /r^) 2 so that the 
emitted energy in the form of the blackbody radiation 
from the heated PC surface may be consistent with the 
energy deposited by the bombardment of gap-accelerated 
e-'s. 

2.7. Effective area of heated polar-cap 

The number flux of e~'s falling onto the polar cap sur- 
face (with area A mp ) can be evaluated by 

f2P* 



N c = k 



2we Agap 



f5.5 x 10 29 kP -2 



M30 



0.2v4. O pc 



(20) 



where A. 



open 



7rr^/t(7LC denotes the entire PC area. 
Since typically 20 % of the open flux tubes thread the 
gap, A gap ~ 0.2A opon is a good compromise. See figure [3] 
for an example of active field lines on the PC surface, 
which indicates that about 20% of the open field line 
fluxes are active. 

From numerical analysis, k ~ 0.3 is appropriate for 
a ~ 75°, while n ~ 0.7 for a ~ 45°. In Newtonian 
approximation, k typically becomes cos a, because we 
obtain B z /B — cos a at the PC surface, where B z refers 
to the magnetic field component projected along the ro- 
tation axis. Thus, k decreases with increasing a. How- 
ever, as the NS rapidly rotates, the space-time dragging 
effect (Muslimov & Tsygan 1992) leads to a smaller k 
than cos a. On the other hand, in a three-dimensional 
pulsar magnetosphere, inward 7-rays preferentially prop- 
agates towards the leading side due to the aberration 
of the photon propagation direction; thus, created cur- 
rent becomes super Goldreich- Julian in the leading side 
to exert a space-charge-limited flow of ions from the 
PC surface in an OG (H06a). Thus, k can exceed the 
value pcj I (OB /27re) at the PC surface. For a highly 
inclined case, the fraction of the field lines having super- 
Goldreich- Julian current density increases because of the 
decreased pgj at the PC surface. Thus, n ~ 0.3 is more 
appropriate than cos(75°) = 0.26 for a — 75°. However, 
for a = 45°, most field lines have sub- Goldreich- Julian 
current; thus, we obtain k m cos(45°) = 0.70. As an in- 
termediate case, we present the solution of k in a three- 
dimensional OG with a — 60° in the left panel of figure[3l 
which shows that the evaluation of k « cos(60°) = 0.5 is 
appropriate on average. 

Each e~ brings the energy ^ra^c 2 on the NS surface, 
where 7h denotes the Lorentz factor at the altitude r = 
rh = v / 27vtolc below which all the emitted photons hit 
the NS. For case 1, most of the e~'s kinetic energy at the 
gap inner boundary will be eventually turned into the 
NS surface emission, because r n /n7Lc ~ 0.028 -^/r^/P 
shows that r nu \\ < r n if P < 0.3 s. For case 2, on the 
contrary, only a portion of jm c c 2 at the null surface will 



be used in the heating of the NS surface, because r nu n > 
Hi- 

For case 1, we can put 7h = 7, where the right-hand 
side is evaluated by equation (|T0|) . Thus, combining 
equations (fT0|) and (|20|) . we obtain 



lh m c c'N c « 4.5 x I0 2 y/ 4 P- 5 / 2 /^ 4 



c 2 l 

1/4 K 



A, 



0.3 0.2A opo 



erg s 



(21) 



For case 2, on the other hand, electrons lose energy by 
curvature radiation from r = r nu n down to r — r n before 
illuminating the NS. Integrating the equation of Lorentz 
factor evolution, 



~dt' 



2e 2 
3P? 



from r = r nu n to r = r n , we obtain 



7h 



2.5 x 10 7 P 1/3 



In 



-1/3 



(22) 



(23) 



where the weak dependence on r nu n/r n , which is of the 
order of unity, may be neglected. Combining equa- 
tions (|2"3")) and l|2"0")) . we obtain 



7 h m e c 2 iV c w7.8 x io 30 p- 5 / 3 



M30 



K A 



0.7 0.2A, 



g a P — 1 

erg s . 



(24) 



Equating the heating rate 7h m e c 2 N e derived just above 
with the emission rate irR pc 2 o~T 4 , where a denotes the 
Stefan-Boltzmann constant, we can evaluate the effective 
area of the heated PC region, ttR 2 c , and obtain 



^ ) = 1.4 x 10- 3 p- 5 / 2 M 30 5/4 
r* 



x 



kT 
100 eV 

for case 1, and 

2 



V 



1/4-7/4 



K A„ 



ap 



* 6 



0.3 0.2A opc 



(25) 



R 



pc 



2.4 x 10" 2 P- 5 / 3 



/ kT 



\ 100 eV 



-4, 



M30 



ap 



'0.7 0.2A opc 



(26) 



for case 2. Note that ^4 ga p/^4opon merely parametrizes 
the fraction of active field lines to the entire open field 
lines. 

2.8. Multiplicity constraint 

For case 1, substituting equation (|2~51 into (ITS|) . we 
obtain 



N y r sa 6.3 x 10" 



n-i 



A v 



0.3 0.2A O p On V 0.67r nuU 



kT 
100 eV 
/ 



(27) 



Death Line of Gamma-ray Pulsars 



5 



For case 2, substituting equation (|2"61 into (fTTj)) . we 
obtain 



x p-35/12 (( 5/4_K Ag a p 



M30 



0.7 0.2A opcn V0.47r nun 



/ 



(28) 



In this paper, we evaluate the NS magnetic moment 
assuming the vacuum dipole radiation formula, /1 2 = 
3Ic 3 PP/(8tt 2 ). Then, equations and give 



-2 / «7 



xP -ll/4p_ 3/4 ^gap_ 



0.3 0.2A opcn 



/ 



(29) 



A 7 r w 1.2 x 10" 1 



kT 
100 eV 



r 5/8 



x p-55/24p_ i5 5/8_ K ^gap 



0.7 0.2A opcn V0.47r nu ii 



(30) 



for cases 1 and 2, respectively. 

Setting 7V 7 r > 1, we obtain the death line on the P 
versus P plane, which describes the condition for an OG 
to emit 7-rays efficiently. Thus, the minimum spin-down 
rate for a given P is obtained as 



lgP= -13.49 + 3.67 IgP, 
lgP= -13.54 + 3.671gP 



(31) 
(32) 



for cases 1 and 2, respectively. 

Since we consider that the gap will not extend in the 
entire open field line region, the obtained death lines (|31|) 
and Q32p gives a tighter constraint than previous works. 
For example, Zhang et al. (2004) gives 

lgP = -14.60 + 3.33 IgP (33) 

for a uniform distribution of a's, and 

lgP = -14.20 + 3.331gP (34) 



for a cosine distribution of a. Death Lines (133)) and (|34l) 
are depicted by the thick and thin dash-dotted lines in 
figure El 

3. DEATH LINE DERIVED BY PAIR-PRODUCTION 
THRESHOLD 

We now turn to the second and important constraint 
obtained by pair-production threshold. To derive this, 
we must evaluate the maximally possible energies of both 
the surface X-rays and the curvature 7-rays. 

First, let us examine the maximum surface temper- 
ature, fcT max , which is probably realized in a limited 
area of the heated PC, near the footpoints of active 
magnetic field lines. That is, we release the assump- 
tion of uniform kT on the heated PC surface, which 
was adopted in § [5J Electrons will fall onto the NS 
surface with the typical Lorentz factor given by equa- 
tion (1231) . where may be replaced with r*. However, 
the logarithmic factor little changes the right-hand side; 



thus, we evaluate the Lorentz factor of in-falling e 's by 
77 = 2.5P 1 / 3 , where 77 = 7/IO 7 . Equating crT max 4 with 
'ym e c 2 N e /A gap = 7m c 2 Ki7P*/(27re), we obtain 

fcT max = 426eV K 1 / 4 7^ /4 p- 1 / 4 M3^n, 6 - 3 / 4 . (35) 
Therefore, the typical X-ray energy can be estimated by 
hv x ~ 3fcT max 

« 0.0023m e c 2 (JL) V4 P- 1 /V 3 o 1/4 n, 6 - 3/4 (36) 
for both case 1 and case 2, apart from the difference of 

K. 

Second, we examine the maximum energy of curvature 
photons. For case 1, we use equation (fT3| with rj — 
2.4 x 10 3 P(r/0.5ro LC ) to obtain 

h Vl = 1.6 x 10 2 m o c 2 p- 7 / 4 ^ 4 r >ti6 1 / 4 . (37) 

For case 2, we compute hv 1 by equation (Ti~4"]) . 

Third, we impose the pair-production threshold condi- 
tion, hux x hv~f > (m e c 2 ) 2 to obtain 



lgP= -14.15 + 2.83 IgP, 
lgP= -14.06 + 2.83 IgP 



(38) 
(39) 



for cases 1 and 2, respectively. 

On these grounds, the actual death line is obtained by 
taking the greater P in each case. For case 1, we obtain 

lgP = max(-13.49+3.671gP, -14.15+2.83 IgP), (40) 

and for case 2, we obtain 

IgP = max(-13. 54+3. 671gP, -14.06 + 2.83 IgP). (41) 

It follows that the death line little depends on a. In fig- 
ure [2 the death line (l40l) is plotted as the (red) thick 
solid line. Case 2 (eq. [H]) is not plotted, because it 
almost overlaps on case 1 (eq. [40]). The (blue) squares 
denote the 7-ray detected pulsars, the (blue) triangles 
the millisecond 7-ray pulsars, and the (blue) circles other 
radio- loud 7-ray pulsars detected with Fermi and/or AG- 
ILE (Aliu et al. 2008; Halpern et al. 2008; Pellizzoni et 
al. 2009; Abdo et al. 2010a; Abdo et al. 2010b; Abdo 
et al. 2010c; Pilia et al. 2010; Saz Parkinson et al. 2010; 
Pilia & Pellizzoni 2011; Saz Parkinson et al. 2011; Guille- 
mot et al. 2011; Keith et al. 2011). The (green) small 
open circles denotes the X-ray detected pulsars (see 
Grindlay and Bogdanov 2009, Becker & Trimmer 1999, 
Becker 2009, for reviews; see also Grindlay et al. 2002, 
Bogdanov et al. 2006, Bogdanov et al. 2010, Bogdanov 
et al. 2011, Saito et al. 1997, Becker et al. 2003, Bassa 
et al. 2004, D'Amico et al. 2002, Eisner et al. 2008, and 
Heinke et al. 2006 for the millisecond pulsars in globular 
clusters; Kuiper et al. 2000, Mineo et al. 2000, Kuiper et 
al. 2002 for J0218+4232; Zavlin & Pavlov 1998, Zavlin 
et al. 2002 for J0437-4715; Bogdanov & Grindlay 2009 
for J0030+0451; Zavlin 2006 for J0437-4715, J2124-3358, 
J1024-0719, J0034-0534 ). The small black dots denote 
radio pulsars that do not show either 7-ray or X-ray pul- 
sations (ATNF pulsar catalog 3 ). 

3 http: / / www.atnf.csiro.au / research / pulsar / psreat / 
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period [s] period [si 

Fig. 2. — Outer-gap death line on the period versus period- 
derivative plane. The (red) solid line denotes the death line ob- 
tained in the present work. For comparison, previous outer-gap 
death lines are depicted as the thick and thin dash-dotted lines 
(Zhang et al. 2004). The thin dotted lines expresses the parameter 
region that gives the labeled magnetic-field strength, -Blo a * * ne 
light cylinder. The (blue) large points denote 7-ray detected pul- 
sars, whereas the (green) small open circles denote X-ray detected 
pulsars (see text for details and references). Left: entire parameter 
space. Right: close up of millisecond parameter range. 

Below the death line, or equivalently in the so-called 
'death valley', the OG activity ceases. It follows that all 
the 7-ray pulsars detected so far and 84 (of 89) X-ray 
pulsars are located above the death line. The remaining 
five X-ray pulsars in the death valley, all of which are old 
nearby pulsars, will be discussed in § 14.21 

It also follows from figure [2] that the new constraint 
obtained by the pair-production threshold (§[3|) becomes 
important when we consider fast rotators (P < 0.1 s). 
To see this, the close-up figure of millisecond parame- 
ter range is presented as the right panel. It shows that 
all the X-ray and 7-ray millisecond pulsars (MSPs) are 
above the death line, including the marginal two cases 
(J1911-6000C in NGC 6752 and J0024-7204C in 47-Tuc). 
If the X-ray emission of such a marginal case is turned 
out to be thermal, it is probably due to the heated po- 
lar cap emission associated with OG activities. On the 
other hand, if they are non-thermal (as in the case of old, 
nearby, non-recycled pulsars), it may be due to a mag- 
netospheric emission from the extended polar gap (see 

4. DISCUSSION 

It is found that the death line has a weak dependence 
on a. This is because negative feedback effects work (Hi- 
rotani 2006) in an OG and its electrodynamic structure 
little changes when we moderately change such param- 
eters as a, P, and [i. However, the gap electrodynamic 
structure is mildly subject to change when we adopt dif- 
ferent kT (eqs. [29] & [30]), because it directly affects the 
efficiency of 7-7 pair production. 

In their death line calculation, Zhang et al. (2004) did 
consider the 7-7 pair-production threshold and wrote 
down the fractional gap thickness, /, as a function of 
a, P, P, and r. In this paper, instead of using /, we 
considered two conditions (-/V 7 r > 1 and pair-production 
threshold) to find that the latter condition gives a tighter 
constraint for MSPs. The difference becomes particularly 
important when we investigate potential MSPs in glob- 
ular clusters, which provide an efficient way to observe 
less-luminous MSPs with long exposures. 



4.1. Death line of millisecond pulsars 

For slowly rotating pulsars (P > 0.1 s), the OG is 
active if the magnetic field strength at the light cylin- 
der, -BlCj is greater than 10 4 G. However, for MSPs 
(P < 0.01 s), the OG has enough luminosity only when 
£?lc > 10 5 G. The reasons are twofold: lower PC temper- 
ature and soft curvature spectrum. For the PC temper- 
ature, equation ([35]) gives fcT max cx BLC 6pll/16 > where 
B LC oc pV2p-5/2 and ^ ^ p i/2pi/2 are ugecL There- 
fore, for the same Plc, the heated PC temperature de- 
creases with decreasing P. For the curvature photon en- 
ergy, equations ^ and (J7J) give hv c oc 

p3/4 p5 /4 There _ 

fore, the maximally possible curvature photon energy, 
hv c , also decreases with decreasing P. 

In spite of the facts above, we should notice here that 
the gap meridional thickness is self-regulated so that the 
pair multiplicity within the gap may be kept unity on av- 
erage. Therefore, the curvature photon energy (or equiv- 
alently, the exponential-cutoff energy) little changes as 
the pulsar ages. This argument is valid until the pul- 
sar approaches the death line, and will be discussed in a 
separate paper. 

4.2. Old nearby pulsars in the outer-gap death valley 

In figure[2] we plot the X-ray detected pulsars by small 
(green) open circles. It follows that several pulsars with 
longer periods (P > 0.5 s) are located in the 'death val- 
ley', which indicates that their OGs have already finished 
activity. They are B0823+26, B1133+16, J0108-1431, 
B0943+10, B0628-28, all of which are old non-recycled 
pulsars with characteristic ages around 10 7 years (Becker 
et al. 2004; Becker et al. 2005; Becker et al. 2006). Al- 
though the photon statistics are limited, B0823+26 and 
B0628-28 are found to show predominantly non-thermal 
spectra. 

We consider that their non-thermal X-rays are emitted 
from the PC accelerator (or polar gap) extended into the 
higher altitudes along the magnetic field lines curving to- 
ward the rotation axis (Fawley et al. 1977; Scharlemann 
et al. 1978; Arons & Scharlemann 1979; Arons 1983). 
For a moderately inclined rotator, such 'toward curva- 
ture' filed lines exist on the opposite side of the OG on 
the PC surface (e.g., in the magnetic azimuthal angle 
range between 105° and 220° in fig. [3] for a = 60°, mea- 
sured counter-clockwise from the +x direction). In an- 
other word, after the OG in the first and fourth quadrant 
ceases activity, the extended polar gap in the second and 
the third quadrant remains active to emit non-thermal 
X-rays or soft 7-rays with small luminosity. Since it is 
reasonable to assume a free emission of e~'s from the 
NS surface by E« (Arons 1981), there will be no 'death 
lines' for this kind of extended polar gap. The space- 
charge-limited flow of e~'s will be continuously acceler- 
ated outwards by a week Eu to emit curvature photons 
from X-ray to soft 7-ray energies, depending on the po- 
tential drop in this gap. The photons will be emitted 
into the direction 9 < a and tt — 9 < a, where 9 denotes 
the photon propagation angle with respect to the rota- 
tion axis. For instance, the total emission solid angle 
becomes 2ir if a — 60°. This corresponds to a mod- 
ification of the original pair-starved polar-cap (PSPC) 
model attempted to apply to the 'away curvature' field 
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lines (Muslimov & Harding 2004a,b; Venter et al. 2009). 
This topic, a modified version of the PSPC model, will 
be discussed in a separate paper. 

4.3. Heated polar cap of millisecond pulsars 

Let us finally consider the relationship with X-ray 
observations. To this aim, we present the OG solu- 
tion obtained for a typical millisecond pulsar parameter, 
P = 10 ms, P — 10~ 19 ss _1 , by applying the numerical 
technique described in Hirotani (20ff). Since the solu- 
tion close to the death line cannot be easily obtained by 
numerical analysis, we consider this set of (P, P), which 
is located relatively away from the death line (fig. [2]). 
To consider a moderate case, we assume a = 60°. We 
present the solved k, the created current density normal- 
ized by the typical Goldreich- Julian (GJ) value, ttB/2ir 
in the left panel of figure El It follows that the created 
current becomes super-GJ, k > cos a = 0.5, in the lead- 
ing side (in the first quadrant), whereas sub-GJ in the 
trailing side (in the fourth quadrant). This is because the 
inward 7-rays propagate towards the rotational direction 
due to the aberration of the photon propagation direction 
to preferentially materialize as pairs in the leading side. 
The created current concentrate near the last-open field 
line in the leading side, because the magnetic flux sur- 
face (i.e., a constant A v surface) has a negative extrinsic 
curvature (like the saddle or the inner surface of a daugh- 
nut) in the lower altitudes where pair production mainly 
takes place. The inward 7-rays preferentially propagate 
towards the lower magnetic colatitudes (i.e., towards the 
magnetic equator) to materialize as pairs near the last- 
open held line. This forms a striking contrast with two- 
dimensional calculation, which predicts that the created 
current peaks in the middle or higher magnetic colati- 
tudes (i.e., away from the last-open field lines) in the 
gap (fig. I of Takata et al. 2008). 

Using k and the Lorentz factors of the in- falling e~'s 
at the PC surface, we can compute the maximum at- 
tainable surface temperature using equation (|35[) . The 
result is presented in the right panel of figure [3] The in- 
tegrated PC X-ray luminosity becomes 3.4 x 10 30 ergs -1 , 
whereas the magnetospheric 7-ray luminosity becomes 
4.4 x 10 32 ergs _1 . We should notice here that this fig- 
ure does not represent the actual temperature distribu- 
tion on the PC surface of a millisecond pulsar by any 
means. For a weak surface magnetic field (B < I0 9 G) 
of millisecond pulsars, the cyclotron energy is compa- 
rable or less than the Coulomb energy. Thus, we can 
expect more or less isotropic heat conduction on the sur- 
face of a millisecond pulsar. The high heat conduction of 
subphotospheric layers, where the in-falling plasmas' en- 
ergy is released, allows the heat to propagate across the 
NS surface, resulting in a greater hot region than given 
in figure [3] (right) with decreasing temperature toward 
the rim of the heated region. Since there is no reliable 
calculation of the heat conduction on a NS surface, we 
simply present the maximum attainable temperature in 
this paper. 

From the X-ray observations of several nearby field 
millisecond pulsars (see Pavlov & Zavlin 1997, Zavlin 
& Pavlov 1998, Bogdanov et al. 2007 for J0437-4715; 
Becker et al. 2000, Becker & Achenbach 2002, for 
J0030+0451; Zavlin 2007 for J2124-3358; Zavlin 2006 for 
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Fig. 3. — Solution of a self-consistent three-dimensional outer- 
gap for an example millisecond pulsar with P = 10 ms, P = 
10 — 19 ss~ 1 , and a = 60°. Both the abscissa and the ordinate 
are in radian unit. The white curve denotes the foot points of the 
last-open magnetic field lines on the polar-cap (PC) surface. The 
magnetic field is assumed to be dip olar centered at the NS cen- 
ter. Left: Distribution of re (eq. |20> . or equivalently, the created 
current density normalized by the typical Goldreich-Julian value, 
QB/(2tt) on the polar-cap surface. Right: Distribution of the max- 
imum attainable temperature [eV] owing to the bombardment of 
the gap-accelerated electrons on the PC surface. Heat conduction 
in the trans-field direction is not considered. 



J1024-0719), it is suggested that the heated PC emission 
can be fitted by a superposition of two thermal compo- 
nents. One component has a higher temperature around 
a few hundred eV with emission area of ~ 0.3 km 2 , which 
is sometimes referred to as the 'core' region of the heated 
PC. Another component has a lower temperature around 
80 eV with emission area of ~ 30 km 2 , which is sometimes 
referred to as the 'rim' region of the heated PC. 

We expect that the higher temperature region appear- 
ing in the leading side (i.e., in the first quadrant of fig- 
ure [3]) near the last-open field line corresponds to the ob- 
served core region, whereas the lower temperature region 
appearing in the trailing side (i.e., in the fourth quad- 
rant) partially corresponds to the rim region. Although 
a detailed calculation of the heat conduction across the 
NS surface may not be easy, it is suggestive that the 
leading side tends to have a higher temperature than the 
trailing side. If the photon statistics will be improved in 
the future, the phase-resolved spectrum of a broad peak 
in the X-ray light curve will show a decreasing tempera- 
ture with pulsar phase. 

This work is supported by the Theoretical Institute 
for Advanced Research in Astrophysics (TIARA) oper- 
ated under Academia Sinica and the National Science 
Council Excellence Projects program in Taiwan admin- 
istered through grant number NSC 98-2752-M-007-006- 
PAE, and the Formosa Program between National Sci- 
ence Council in Taiwan and Consejo Superior de Investi- 
gaciones Cientificas in Spain administered through grant 
number NSC100-2923-M-007-001-MY3. 

5. APPENDIX 

Derivation of equation (j4|) is given below. For simplic- 
ity, we assume a Newtonian dipole magnetic field config- 
uration. We consider the plane (x,z) in which both the 
rotation and magnetic axes reside. Let us express the 
position of the cross section (point PI in fig.[T]) between 
the null surface and the last-open field line as (xi,zi) in 
Cartesian coordinates. Then, the radial distance of point 
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PI from the NS center becomes 



n - a) 



sin (6»lc - a) 



(42) 



where nu n specifies the null surface, and #lc the co- 
latitudes (measured from the rotation axis) of the point 
at which the last-open field line becomes tangent to the 
light cylinder. Thus, the distance of point PI from the 
rotation axis becomes 



For photons emitted from the middle colatitude in the 
gap, we thus obtain the propagation distance below the 
critical field line, 



I ~ Xl ^ 2 = 7^ ( sin ^ nul1 ~ tan 02 cos 0nu11 ) ( 46 ) 



sm (tf nul i - a) . 



sin (0 LC - a) 



(43) 



where 

cos a tan 9 2 — sin a 
1 + tan 2 6-2 



cos( 



all sm 



7 LC 



a 



= e. (47) 



If a photon is emitted leftward from point PI, it crosses 
the critical field line (fig. [T|) before escaping from the gap 
at point P2 (x 2 ,z 2 ). For a Newtonian dipole field, we 
obtain 



r-2 



? ■> sin (0 2 - a) 
.x 2 + z 2 = n7 LC — -2-^ ^ , (44) 



J (^null - a) 



and hence x 2 = r 2 sm9 2 . Therefore, we obtain the dis- 
tance between PI and P2, 



xi — x 2 — ri(sin( 



null 



tan^ 2 cos^nuii) (45) 



The factor \[2 in the denominator of equation (|4"5)l is 
appropriate if the poloidal magnetic field lines are co- 
centric. For a dipole field geometry on the poloidal 
plane, equation (|4"5f overestimates I. Nevertheless, the 
inward 7-rays propagate towards the leading side in a 
three-dimensional magnetosphere owing to aberration 
and propagate a longer distance than what is evaluated 
only on the poloidal plane. Thus, equation ([4T>| gives 
a good estimate. Since e <C 1, tan #2 ~ tana reduces 
equation (14T>1) into equation 
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